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We study a simple model for the metallic stripes found in Lais-^NdoASr^CuOi: two chain Hub- 
bard ladder embedded in a static antiferromagnetic environments. We consider two cases: a "topo- 
logical stripe" , for which the phase of the Neel order parameter shifts by tt across the ladder, and a 
"non-topological stripe", for which there is no phase shift across the ladder. We perform one- loop 
renormalization group calculations to determine the low energy properties. We compare the results 
with those of the isolated ladder and show that for small doping superconductivity is enhanced in 
the topological stripe, and suppressed in the non-topological one. In the topological stripe, the 
superconducting order parameter is a mixture of a spin singlet component with zero momentum 
and a spin triplet component with momentum tt. We argue that this mixture is generic, and is due 
to the presence of a new term in the quantum Ginzburg-Landau action. Some consequences of this 
mixing are discussed. 



Recent neutron experiments have demonstrated that 
for Lai Q^^NdoASr^CuO^, doped holes segregate into 
an array of "stripes" embedded in an antiferromag- 
netic background.u Moreover, it was found that for x — 
0.12,0.15,0.20 the systemspexhibit simultaneous stripe 
and superconducting order .□ According to Tranquada et 
al, a consistent interpretation of the observed spin and 
charge incommensurability is that the stripes are the an- 
tiphase domain walls of the magnetic order. A number of 
theoretical works have addressed issues varying from the 
origin of the stripe oijdprOl, to the effects of stripe order 
on superconductivityu'u. 

In this paper we examine a simple model of a single 
stripe. The model consists of a Hubbard ladder (where 
the holes reside) embedded in one ofjlhe two types of 
static antiferromagnetic environments. □ In case (a), the 
magnetic order on the two sides of the ladder forces a tt 
phase shift in the Neel order parameter (Fig. 1(a)). We 
shall subsequently refer to this type of stripe as being 
topological. In case (b), the Neel order parameter stays 
in-phase across the stripe (Fig. 1(b)). We refer to this 
case as being non-topological. 

The low energy properties of the model are determined 
by performing one-loop renormalization group (RG) cal- 
culations. Comparing our results with the results of sim- 
ilar calculations for an isolated ladder allows us to infer 
the effects of a magnetic environment on the low energy 
properties of the two chain Hubbard ladder. We empha- 
size that in our view the present work sheds light on how 
superconductivity survive stripe ordering, but not how 
stripe ordering trigger superconductivity. 

Our main results are summarized as follows. When 
the Hubbard ladder is lightly doped, we find that super- 
conductivity is enhanced (relative to the isolated Hub- 
bard ladder) in the topological stripe, and is totally sup- 
pressed in the non-topological one. The superconducting 
order parameter in the former is a linear combination 



of a spin singlet component with zero-momentum, and 
a spin triplet component with momentum tt. We argue 
that this mixture is generic, and is due to the presence of 
a new term (Eq.(p^) in the quantum Ginzburg-Landau 
action. 
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FIG. 1. The topological (a), and non-topological (b) 
stripes. The arrows represent the magnetic moments of the 
environment. Antiferromagnetic coupling exists between an 
arrow and electrons on the site next to it. To determine 
whether a stripe is topological, one needs to interpolate the 
magnetic order across the ladder. 



The model 

The model we shall study describes a metallic "river" 
(stripe) embedded in an insulating antiferromagnetic 
background. The stripe will be modeled by a "Hubbard 
ladder" (Fig. 1). The couphng between the electrons in 
the stripe and the magnetic moment of the background 
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is the usual antiferromagnetic spin-spin interaction. The 
Hamihonian is 

Ha = -tJ2{[4+i^('"^'T + c^^Kl^c^^^ + c+d;^ + h.c] 
-m[c+ + dtM + Ma(-l)'[c+ + 7ld+d,„]}, 
^ = ^ YlytaC^''^^-aC^-'y + d+d,„d+_^di^„]. (1) 

In the above, i runs through the lattice sites, a = ±1 
labels the electron spin, c and d annihilate electrons on 
the two chains of the ladder, /i is the chemical poten- 
tial, U is the Hubbard interaction, M is the internal 
staggered magnetic field induced by the background mo- 
ments, and 77 = -|-I(— 1) when the stripe is topological 
(non-topological) . 

When M = Eq.(|l|) describes an isolated Hubbard 
ladder. The following is a brief summary of the theoret- 
ical results for this system.El At half filling, the ladder is 
a Mott insulator with both a spin and charge gap. Su- 
perconductivity develops when it is lightly doped. In the 
superconducting phase, the order parameter is the out-of- 
phase linear combination of the usual spin singlet order 
parameters in the symmetric and antisymmetric bands. 

When M ^ Q the first Brillouin zone is halved, and 
the free-particle bands of the isolated ladder hybridize. 
The resulting bandstructures are shown in Fig. 2. At half 
filling there are four Fermi points (which pairwisely co- 
incide) for T] = +1 (Fig. 2a). In contrast, an energy gap 
A ~ M opens up for 77 = —1 (Fig. 2b). Doping moves 
the chemical potential so that for both 77 = ±1 there are 
four distinct Fermi points (fc^o ± Skp)- Here Skp is the 
shift in Fermi momentum due to doping, and kpo satis- 
fies {M/t)^ +Acos^{kFo) = 1 for ry = 1, and kpo = ±7^/3 
for 77 = —1 respectively. (In the above and hereafter we 
will set the lattice constant to unity.) The dispersion of 
the energy band that intersects jj, at, e.g., kpQ + Skp is 
such that dEk/dk < at the Fermi level. We shall define 
kpfi- = kpQ + Skp, and refer to the portion of energy 
band in its vicinity as the (i?— ) branch. (Here R desig- 
nates the right group of Fermi points, and — reflects the 
fact that dEk/dk < 0.) Similar notation will be given to 
the Fermi momenta and energy band branches associated 
with the other three Fermi points. 

Let ipR± and ipL± be the electron annihilation oper- 
ators associated with (-R±) and (L±) branches respec- 
tively. In terms of them the lattice annihilation operator 
Cio- is given by 

= e*'=-«-^'7/.ii_.(x,) + e**'--«+^'Vi?+.(^») 

+ e*'=---"'VL-.(x.) + e*^---+^' Vl+.(xO. (2) 

In the above Xi is the coordinate of the i'th lattice site. 
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FIG. 2. The bandstructure of the topological(a) and 
non-topological(b) stripe, fio is the chemical potential at half 
filling; /ii is the chemical potential after doping. 

The RG solution 

Since we are primarily interested in the low energy 
properties, we shall limit ourselves to a small window of 
the energy band near each Fermi point. To be more pre- 
cise, we assume that we can truncate the Hilbert space 
so that aside from the free Fermi sea, it only includes 
the subset of particle-hole excitations where the parti- 
cle/hole states fall inside the small windows. Further- 
more, within each window we shall ignore the curvature 
of the band. Next we project the Hubbard interaction 
onto the truncated Hilbert space. The result is a family 
of vertex functions which describe the scattering matrix 
elements between states near Fermi energy. In order to 
simplify the calculation, we shall assume that the ver- 
tex function has -f <-^- — symmetry. This requires that 
Vp^ =1 dE^^ /dk I is the same for each energy band E^\ 
In the rest of the paper we shall concentrate on small 
doping where this symmetry is approximately valid. 

For generic doping, i.e. when there is no umklapp 
scattering, there are five independent vertex functions 
given by 

+54aC^lC^2(5<T 
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All other vertices that can be generated from the above 
r's by letting L and/or + ^ — have the same value. 
These vertex functions are characterized by 18 coupling 
constants gl^ and gl^. In the above we have chosen the 
notation so that i and j reflect the changes in the upper 
(+ or — ) and lower (i? or L) indices during scattering. 
The convention is that 1 =back scattering, 2 = inter- 
branch forward scattering, 3 — umklapp scattering, and 
4 = intra-branch forward scattering. Moreover, s and a 
in the subscript of g charactcrizt; the spin dependence of 
the scattering matrix elements. □ 

The goal of the RG calculation below is to progres- 
sively eliminate high-energy electronic excitations, and 
study its effects on the remaining lower-energy ones. In 
this work the mode elimination is implemented pertur- 
batively. To be more precise, we integrate out those 
'>Pijcr{kys [i = L/R,j = +/—) where IL — dE<VF \ 
k \< Ec within one- loop approximation.^ (Here i?c is a 
microscopic energy cut off.) Since the calculation is per- 
turbative in nature, the validity of the results relies on the 
smallness of the g 's. Straightforward calculations show 
that among the 18 coupling constants only 12 renormal- 
ize. The resulting recursion relations are given by 
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In the above dZ = ^f^liil"^)? where i? is the running 
energy cutoff. In order to apply these recursion relations 
we need to determine the initial values of g's. It turns out 
that the latter strongly depend on the topological type 
of the stripe. Let us concentrate on the case of small 
doping where kFR+ ~ kpR- and /c_fl+ ~ kpL-- 
For the topological stripe ( = 1 ) we have 
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g^a = gia = 52a = 4a-5 U, 

gla = gla = gla = o> 
gl 



gl 
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For the non-topological stripe ( 77 = — 1 ) we have instead 



2 2 2 rr 

54s = 5ls = 52s = U, 



'^gla = gla = '^gla = 52a = gla = "^gla = ^o'h'U, 

g\s = 52s = gls = (1 - la'b^W. (6) 
In Eqs.d) and (§) 
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In Eq.(||) A = 2cos(fcFo) for 77 1, and A = 1 - 
2cos{kFo) for JI= —1 respectively.t2l 

Given Eqs.(^ and (^), we numerically iterate Eq.(Q) to 
determine the renormalized (/'s. The result is trustworthy 
when i) U/t << 1 and ii) all renormalized g's are << 1. 
Under i) and ii) we find that in all cases there exist more 
than one g that grow upon renormalization. To deduce 
the implication of these growing coupling, we compute 
physical susceptibilities. 

For [7 = 0, the following 12 susceptibilities capture 
all logarithmically divergent ones up to i? <-!■ L and/or 
-|- <-> — exchanges: 



Xa(^) = J dxdte''^' < T[0„(x,t)0+ (0,0)] >, (8) 
where 

Ocdwiix) = i(;lf^^{x)i(;R^'^{x) -I- '0|j^|(a;)V'i?.-i(a;), 
Osdwi\\{x) = [ip'l^^{x)ipR^^{x) - 'iljl^^^{x)ijjR^i{x)]/2, 
Osdwi±{x) = [ijj]^^-^{x)i/jR^i{x) +-il;l^^^{x)ipR^^{x)]/2, 
Osslix) = [?/'Ji+|(a;)V'Jj_|(a;) - V']^,+x(a;)^l?._|(a;)]/\/2, 
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Otsi±{x) = 2pl^^^{x)ipl^_^{x), 
Ocdw2{x) = ■iIj]^^^{x)^/jl-t{x) + Vjj.+^(x)V'L_x(a;), 
Osdw2\\{x) = [tpli^^{x)ipL-^T{x) - V'|j^|(a;)V'L-i(2^)]/2, 
Osdw2±{x) = [-0]j^^(a;)-0L_x(x) + tpl^_^_^{x)tpL--\{x)]/2, 
Oss2ix) = [i;]^^^ix)i:[_^{x) - Vj]^,+|(a;)Vj[_^(a;)]/\/2, 
Ots2\\{x) = [ip]i+^{x)i;l_^{x)+ij]^^^{x)ijl_^{x)]/V2, 

Ots2±(x) = 1p'l^-^{x)lpl^^{x). 

(9) 

In the above sdw, sdw label the charge and spin density 
wave susceptibilities, while ss and ts label the supercon- 
ducting ones. To the lowest order in g's, the results are 

XaM =^aX0M( — XQ = ^H — )- (10) 

In the above the amplitudes A are given by Acdwi = 
Acdw2 — 2 and Aa ~ 1 for all other a's, and the expo- 
nents Xa are given by 

— 2 I 2 _ 9 1 
Xsdwiw = gls + .94a - 2g4„, 

2 2 
a^sduili = 54s ~ 54aJ 

_ 2 1 I 2 I 1 

2^ssl — —54s — 54s + 54a + 54a J 

_ 2,1,2 1 
Xtsl\\ — —9is I .94s ' 54a — 54a J 

xtsii^ = -54s + als - gla + ala^ 

Xcdw2 = 52s - '^als+gla^ 
Xsdw2\\ = gls - 2.9ia + 52a: 

Xsdw2i^ = gls - glai 

Xss2 = -gls - g\s + .92a + 5^, 
Xts2\\ = .9is - 52s - 5la + 52a > 

a:ts2_L = g\s - gls + gla - gL- (n) 

By substituting the renormalized g's into Eqs.(|o|) and 
( ITll ) we can calculate Xa{^) as w ^ 0. In the following we 
shall restrict M/t << 1. We find that, like the isolated 
ladder, Xss2 is the most divergent susceptibility for 77 = 1. 
For 77 — —1, we find that Xcdwi and Xsdwi± are the most 
divergent while Xss2 is non-divergent. 

By setting^^^ = we reproduce the results for the iso- 
lated ladder .0 After a comparison, we find that for 77 = 1 
superconductivity is enhanced by the non-zero g^^. In 
other words, superconductivity is enhanced in the topo- 
logical stripe! 

Before proceeding to the symmetry of the supercon- 
ducting order parameter, we comment on the two condi- 
tions under which the above results are obtained, namely, 
M/t << 1 and small doping, i) The dependence on 
M/t: for 77 = 1 we find that superconductivity is re- 
placed by spin density wave {SDWl _L) when M/t ex- 
ceeds ~ 0.4. On the other hand, for rj — —1, Xcdwi and 
Xsdwi± reamins as the most divergent susceptibilities for 



all values of M/t. ii) The dependence on doping: for 
large doping, the symmetry of -I- <-> — is lost. In that 
case the number of independent vertex function greatly 
increases. We have not done the RG calculation for the 
general case. What we say is that for a fixed M ^ 
the effect of environment diminishes upon doping. Thus 
at large doping the results should resemble that of the 
isolated ladder. For the latter it is found that supercon- 
ductivity survives for doping up to sa 50%. 

Symmetry of the superconducting order parame- 
ter 

As XSS2 diverges the following operators tend to de- 
velop expectation values: 

Kix) = i'R+^ix)lpl^iix) - 'lljl^^^{x)ijl_^{x), 

bti^) = i'l-^ix)i'l+^ix) - ^^j,_^ix)^[+^ix). (12) 

A straightforward mean-field analysisQ indicates that 
the out-of-phase combination of 6^ and 6^, i.e., 

B+ = I dx[bt{x)~b+{x)], (13) 

acquires non-zero expectation. Of course, after including 
phase fluctuation, we expect < [x)B{x') > to have 
only quasi-off-diagonal long-range order even at T = 0. 

bi and 62 might appear to be singlet Cooper pair anni- 
hilation operators. In fact this is not so. Indeed, due to 
the term proportional to M in Eq.(^, spin up and spin 
down electrons experience different one-body potentials. 
Consequently, and V'A-i-i differ by more than a spin 

flip. As a result, the Cooper pair created by i?+ has both 
singlet and triplet characters. To see that more clearly, 
we express ^~^j^^{k) and ip'j^j_^{k) in terms of the band 
operators at M = 0: 

i'R+aik) = ac+„(fcFfi+ + fc) + abc+„{kFR+ + fc - tt), 
i>R-A^) = ac+^ikpR- +k-Tr)- abc+^ikpR- + k), 
i't+ai^) = ac+^{kFL+ + fc + tt) - abc+„{kFL+ + k), 
^t-aik) = act^ikpL- + fc) + abct^ikpL- + fc + tt). (14) 

In the above c^^ and c^^ create an electron in the sym- 
metric and antisymmetric band respectively. Substitut- 
ing the above results into Eq.(p^ we obtain 

/dk 
-{a2[c+(Kfc))c+(-Mfc))-(T-i)] 

-6^[c+ (p(fc) + ^)c+ (-p(fc) - tt) - (t^i)] 
+ab[c%{p{k) + ^)c+ (-p(fc)) + (T^i)] 
-a6[c+ (p(fc))c+ (-p(fc) -n) + (t-i)]}, 

/dh 
-{a2[c+ (g(fc) + 7r)c+ (-q(fc) - ^) - (T^i)] 

-62[c+(9(fc))c+(-g(fc))-(T-i)] 
-a6[c+(g(A;))c+(-g(fc)-7r) + (t^i)] 
+ab[cUq{k) + 7T)ct^{~q{k)) + (T-^i)]}. (15) 
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In the above p{k) = kFR+ + k and q{k) = kpR- + k. In 



Eq.(|15D the terms proportional to a and h have S — 
and momentum = 0, while those proportional to ah have 
5 = 1 and momentum = tt. 

A Ginzburg-Landau theory for the mixed sym- 
metry 

The result obtained above for the mixed component 
order parameter can be understood on the basis of more 
general symmetry considerations. For a system hav- 
ing low energy antiferromagnetic and superconducting 
modes, the following term is allowed in the effective ac- 
tion: 

5™„ = A / dtd'^xN{x,t)-['$t{x,t)^*{x,t)+c.c]. (16) 



In Eg. ([l6|) , N is the Neel order parameter, is a zero- 
momentum spin singlet Cooper pair field, and "0* is a 
momentum Q = (tt, ...,7r), spin triplet Cooper-pair field. 
In addition to satisfying the translational symmetry, the 
fact that the triplet Cooper field has a center-of-mass 
momentum Q is crucial for Eq.(|l^) to be allowed by the 
point group symmetry. Let us assume that 



k 

•' k 



d^xii^t) 
where 



„+ 



fc+QT -fei 



(17) 



0(fc) = cf,{~k), 
xik)^-xi-k-Q). 



(18) 



The fact that N transform as the identity under the ac- 
tion of the point group about a lattice site constrains the 
symmetry of (f>{k) and x{k)- In two dimensions, where 
the point group is Abelian, (/)(fc) and xik) must transform 
identically: 

X(fc)cx0(fc). (19) 

Putting Eqs.([T^) and (^9|) together we obtain 

0(fc) = -0(fc + Q). (20) 

This can be satisfied, e.g., by either for dxi_y2 pairing, 
where 0(fc) oc cos{kx) — cos{ky) , or for anisotropic s-wave: 
(^(fc) cx cos{kx) + cos{ky). 

Depending which one among {N ,ipt,fp} acquires ex- 
pectation value, Smix causes mixing between the remain- 
ing two. For example, if < N >^ 0, Smix mixes "0 and 
■04. This is preciseh2-| what was pointed out by Schrief- 
fer. Wen and Zhangl3 and is a result of our calculations. 



As another example, if < -0 0, then N hybridizes 
with ipf As a consequence, after integrating out ipt, the 
Greens function of N is modified so that 



G„n{q,uj) 



1-A2 I V'PG„„(g,a;)G„(g,c^)■ 



^21) 



In the above G„„ and Gu are the Greens functions asso- 
ciated with N and ^pt respectively. A resonance appears 
when 



1 



(22) 



Recently Demler and ZhangEl postulated that Gtt ex- 
hibits an isolated pole, and that the solution of Eq.(|2^) 
near this pole is the observed resonajice in the neu- 
tron scattering of several high Tc oxides.Ej We emphasize 
that i) the hybridization between N and ■0t predicted by 
Eq.(|l6|) is generic. It does not rely orv-a|n enlarged symme- 
try group for the order parameters.tS ii) The resonance 
in can also appear near the resonance in G„„. In 
that case the mode is predominately magnetic in origin. 
This would explain the facts that the observed cross sec- 
tion and the Cu — Gu bilayer modulation of the resonant 
mode are very close to those of antiferpamagnetic spin 
waves in the undoped antiferromagnetJlj iii) The cou- 
pling described by Eq.(|l6|) exists in the normal phase 
as well. As before, this coupling modifies G„„- Since 
Gsa is no longer of the form | i/j p (5(w)(5(g)|-a.bove T^, it 
can cause the broadening of the resonance.EJ Of course 
in the latter circumstance, quasiparticle excitations will 
also contribute to this brodening. 

Some Concluding Remarks 

According to Tranquada et al, superconducting and 
stripe order coexist for Lai,e-xNdoASrxGuOi with x = 
0.12, 0.15, 0.20. It is tempting to relate a single stripe in 
such system to the model we considered above. There 
are two immediate worries in drawing such a connection. 
The first is that a real stripe is not necessarily made up 
of a pair of chains. The second is that a real stripe will 
not be absolutely straight. 

We do not think the width of the stripe (as long as 
it is reasonably small) will qualitatively change our con- 
clusions. The basis for that belief is that it is known 
for the n-leg Hubbard ladders that away from half filling 
the value of n does not affect the £act that the system is 
superconducting for small doping.li3 

Now we come to the shape of stripes. When the stripes 
meander, the potential seeing by the electrons is no 
longer invariant under the translation by two lattice spac- 
ing parallel to the stripe. The irregularity presents itself 
as both scalar and magnetic impurities. They both have 
pair-breaking effect on the superconductivity as should 
be the case for a "d-wave" order parameter. Nonethe- 
less, the basic fact that local Neel order mixes singlet 
with triplet remains true. 

No comparison with real systems can be made without 
understanding the effects of coupling between stripes. In 
the absence of such coupling, superconducting long-range 
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order is not possible. Knowing how Cooper pairs tunnel 
from one stripe to another is crucial for the understand- 
ing of, e.g., the experimental Tc versus 5q relation (here 
5q is the incommensurability in neutron peak). Presum- 
ably tunneling is most efficient when two stripes "collide" 
with each other. If one assume that the average distance 
(along the stripe) between the collision points scales with 
the average distance between the stripe, then Tc will scale 
linearly with 5q as experimentally observed. _ 

In a recent paper, Emery, Kivelson and Zachaiu at- 
tribute the superconductivity in high compounds to 
the "spin gap proximity effect" . In that mechanism, 
due to the tunneling of a pair of electrons from a sin- 
gle metallic chain to the environment and back, a spin 
gap is induced on the former. In that case, because of 
the spin-charge decoupling in the chain, the charge slid- 
ing mode is left as the only low energy degree of freedom. 
In Ref.u this is identified as superconductivity. The pic- 
ture emerges from our study differs somewhat from the 
above. The biggest differnece is that the environment 
is not the cause of superconductivity. To be more spe- 
cific in our model i) The magnetic environment is not 
spin gapped. ( We note that experimentally the qua- 
sistatic spin peak has been observecUin the "striped" 
compounds in Lais^xN daASrxCuO^) ii) The stripe 
is superconducting even if the coupling to the environ- 
ment is switched off. We should, however, mention that 
the coupling to the enviroment did enhance superconduc- 
tivity. Finally, we point out that in this work we have not 
addressed the issue of long-range Coulomb interaction on 
the superconductivity in stripes. 
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